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Abstract. Recent studies of galactic nuclei suggest that 
most galaxies are seeded by super-massive black holes 
which power the central nucleus. In this picture, the proto- 
galactic object is likely to have undergone a very active 
phase during which the surrounding medium was shocked 
and heated up to very high temperatures. We predict the 
cosmological effects of this scenario in terms of the ther- 
mal and kinetic Sunyaev-Zel'dovich distortions induced on 
galactic scales by a population of proto-galaxies. These 
predictions are compared to the observational limit on the 
mean Compton distortion set by the COBE-FIRAS instru- 
ment. This enables us to derive tight constraints on the 
fraction of proto-galaxies seeded by super-massive black 
holes as well as on the black hole-to-spheroid mass ratio. 
Finally, we estimate the contribution of such a population 
to the angular power spectrum of the Cosmic Microwave 
Background temperature anisotropics on very small angu- 
lar scales (I ~ 10 4 - 10 5 ). 

Key words: Galaxies: formation, Cosmology: cosmic mi- 
crowave background 



1. Introduction 

Most theories of hierarchical structure formation are based 
on the study of the evolution of density perturbations un- 
der their own gravity. A density fluctuation, which repre- 
sents an over- (or under-) density with respect to the mean 
matter distribution, contains both baryonic and dark mat- 
ter (DM). The baryonic component sinks into the gravi- 
tational potential of the DM halo. It collapses and cools, 
resulting in star formation. In these scenarios, after the 
gravitational collapse of the DM halo, stars are assumed 
to be the first objects to form. A structure will thus end 
up as an emitting object after virialisation has occurred. 

An alternative picture involves the formation of 
a super-massive black hole (BH) that pow ers the 
central regions of galaxies Lyndcn-Bcll(1969) , . Nu- 



quasar luminosity function to galaxy formation scenar- 
ios by assuming that the formation of quasars (i.e. 
BH) in the potential well of the DM halos consti- 
tutes one of the phases in the galaxy formation pro- 
cess [Efstathiou k Rees(1988)[ [Haehnelt k Rees(1993) 



Nusser k Silk(1993), Haiman k Loeb(1997), . Recent ob- 



servations even suggest that a super-massive BH may be 
present in the centres of all galaxies with spheroidal com- 
ponents Kormcndy k Richstonc(1995)[ . 

Several authors have looked at several consequences 
of the presence of massive BHs on galaxy formation and 



evolution 


Haiman k Loeb(1997) 




Natarajan ct al.(1998) 


f3ilk k Rccs(1998) 


. In this paper, 


we investigate the cos- 



mological implications of such an alternative scenario for 
the Cosmic Microwave Background (CMB) anisotropics 
and spectral distortions. More specifically, we study the 
effects of the outflows, driven by the BH activity, on the 
gas within the seeded proto-galaxy. In fact, the outflow 
expands and shock-heats the ambient medium (proto- 
galactic gas), and then interacts with the inter-galactic 
medium (IGM). Three regimes of interest may be consid- 
ered: 1) the high density region of the proto-galaxy, 2) the 
low density IGM and 3) the thin compressed layer (four 
times denser than the IGM) induced by the front shock. 
The second and third regimes give results very similar to 



those computed in Aghanim ct al.(1996). We thus focus 
on the first regime, i.e. the localised effects of the BH- 
driven shock on the gas within the seeded proto-galaxy. 
This shock-heated gas will Compton scatter the CMB 
photons and induce spectral distortions and temperature 
anis otropics through the so-calle d Sunyaev-Zel'dovich ef- 
fect [Sunyacv k Zel'dovich(1980)| , . The thermal SZ effect 
depresses the CMB brightness in the Rayleigh- Jeans re- 
gion and increases it above a frequency of about 219 GHz. 
Its amplitude represents the integral along the line of sight 
of the electron pressure. It is proportional to the electron 
density n e and is characterised by the Compton parameter 
y defined by: 



merous studies have been performed that relate the 
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where T is the temperature of the gas, R the physical 
size of the structure, m e the electron mass, k the Boltz- 
mann constant, c the speed of light and ctt the Thomson 
cross section. An additional secondary anisotropy arises 
due to the first-order Doppler effect of the CMB photons 
when they scatter on a structure moving with respect to 
the Hubble flow, with radial peculiar velocity v r . This in- 
teraction is called the SZ kinetic effect. It generates an 
anisotropy, with no specific spectral signature, whose am- 
plitude is given by: 




(2) 



Previous work on galaxy formation has evaluated 
the global distortion of a population of galaxies in the 
virialised regime. The global Compton parameter was 
found to be much smaller than the constraints set by 
the FIRAS instrument (Far InfraRed Absolute Spectro- 
photometer) on board COBE (COsmic Background Ex- 
plorer) Fixsen ct al.(1996), on the global SZ distortion of 
the universe. By contrast, our model focuses on a regime 
in which proto-galaxies undergo a BH formation phase 
that induces larger distortions. The paper is organised as 
follows: in §2, we model the shock in an individual struc- 
ture and give its physical characteristics (size and temper- 
ature). In §3, we compute the predicted number density 



of primordial galaxies, using the Press & Schcchter(1974) 
mass function. In §4, we generalise the description of the 
shock to the whole population of proto-galaxies and we 
simulate maps of the induced secondary anisotropics. We 
estimate this contribution to the CMB anisotropies. We 
also compare our predicted global y parameter to the 
COBE-FIRAS value and derive constraints on the model. 
Conclusions are given in the last section. 

2. Modelling the shock 

In the galaxy formation canonical scenario, a galaxy forms 
in the gravitatio nal potential well o f a DM halo of mass 
M ha i . Following jSiik fc Roes (1998] , we consider that each 
forming galaxy, which is fully described by the mass 
Mhalo, hosts a super-massive BH. The fraction of seeded 
proto-galaxies is thus 100%. For the sake of simplicity, we 
assume that all proto-galaxies are spheroids. A more re- 
fined and accurate description would take into account 
morphological segregation Balland et al.(1998), . How- 



ever, these effects do not introduce significant differences. 
We assume that the BH radiates during its lifetime tsH 
a fraction 6e = 0.1 — 0.2 Natarajan & Sigurdsson(1999), 



of its Eddington luminosity (Lbh = ^E^edd)- The latter 
is fixed by the BH mass, Mbh'- 

AirMBHrripGc 



L e dd = 



(3) 



where m p is the proton mass and G the gravitational con- 
stant. The mass of the BH should be directly related to 



the mass of the proto-galaxy it seeds, and thus to the 
fraction of baryonic matter locked up in the spheroid 
in which the central BH collapses. We assume a sim- 
ple relation of proportionality (Mbh = £bh M sp h) be- 
tween the mass of the BH and the mass of the spheroid 
Hachnclt fc Rces(1993)| give e B H ^ 10~ 3 
that ebh declines with 



M B ph 

they assume 

observations of galactic nuclei Kormendy et al.(1997) 



mass. 



and 
However, 



Magorrian et al.(1998), indicate that cbh is relatively 
constant. Here, we take the value €b h — 2.10~ 3 as ad- 
vocated by [Magorrian ct al. (1998)1 and [Silk fc Recs(1998 



We assume that the medium is instantaneously ionised, 
which is very likely due to the intense radiation emitted 
by the central BH Voit(1996), . We therefore study the 



propagation of a strong shock driven by a mechanical en 
ergy which represents some fraction of the BH luminos 
ity. The first analytic solutions of spherically symmetric 
explosions were given by Taylor (1950) and |Sedov(1959) 
and applied to several astrophysical problems such as stel- 



lar winds and supernovae explosions Ikeuchi et al.(1983) 

Bcrtschinger(1986)|, |Koo fc McKee(1992a) 

Koo fc McKee(1992b)| , |Voit(1996)| , . Following |Voit(1996) 



we study the expansion of the shock and its effects on 
the gas within the seeded proto-galaxy. We characterise 
the shock in an expanding universe by a set of physical 
properties (its radius, velocity and temperature). Adapt- 
ing Voit's solutions to our study, the shock at a redshift 
z < z», where z» is the redshift at which the BH switches 
on. has a radius R s which reads as: 



R s 



ttGE 



1 + VL z 
1 + Qqz* 



1/5 



(1 + ^qZh.) 1 / 5 
(1 + Z*) 2 / 5 



1/2 



-i 2/5 



(1 + Z)~\ 



(4) 



It propagates at a velocity 



nE Q H G ^(l + z*) 3 



3sn b (i + iv*) 3/2 



1/5 



1 + Sl z 
1 + Qqz* 



1/2' 



-i -3/5 



l + Z 
1 + Z* 



(5) 



In the two previous equations, Eq = eLbh^bh represents 
the mechanical energy which drives the shock. Very little 
is known about the fraction e which can be on the order 
of 0.5 or even higher Natarajan & Sigurdsson(1999), . For 



the lifetime of the BH in the bright phase, we choose the 
recent value derived by Haiman & Loeb(1997): Ibh = 10 6 
years (our results, however, are not very sensitive to the 
exact value of tsn)- <5 is the mean over-density of the 
proto-galaxy (assumed to be spherical), H is the Hubble 
constant, tto is the density parameter and fib is the baryon 
density in the universe. Zero subscripts denote present 
day values. In the following and throughout the paper, 



3 



we use 



Walker et al.(1991 



ff /(1 00km/Mpc/s) = 0.5 and Q b = 0.06 
and give the results as a function of 



spheres of 8h 1 Mpc, erg = 0.6 (as cluster-normalised, e.g. 



f2o, for fio = 1 an d = 0.3. 

Following Natarajan fc Sigurdsson(1999)[ we assume 
that the BH radiates 10% of its Eddington luminosity, 
i.e. 6e = 0.1 and that half of the corresponding energy is 
mechanical. Furthermore, we consider that galaxies form 
at sufficiently high redshifts (> 5) so that only hydrogen 
and helium are present in substantial amounts. We com- 
pute both the size and velocity of the shock and find that 
the host proto-galaxy is always embedded in the shocked 
region. 

2.1. Cooling mechanisms 

The temperature of the shocked medium can be directly 
derived from the equipartition of energy: 



(6) 



Here, /i = 0.6 is the mean molecular weight of a 
plasma with primordial abundances. The galactic mat- 
ter can be shock-heated up to very high temperatures. 
For the most massive galaxies, we find that the mat- 
ter is heated up to a few 10 8 K, a value compara- 
ble to the temperature of the intracluster medium in 
galaxy clusters. For these temperatures and redshifts (> 
5), the main cooling process at play is bremsstrahlung. 
Therefore, the shocked gas looses heat with a cooling 
rate given by a temperature-dependent cooling func- 
tion A(T). We have compared three different cool- 



ing functio ns Bertschinger(1986) Koo fc McKee(1992a 



Voit(1996) , given in the literature for temperatures be- 
tween a few 10 5 and 10 8 K. We find that the condition 
required for efficient cooling is always satisfied in our red- 
shift range. This means that the cooling time, given by 
tcooi — 5kT/nn A(T) with n B the hydrogen number den- 
sity, is always smaller than the age of Universe. Cooling is 
thus very efficient in our picture. The comparison between 
the three cooling functions shows that they all give essen- 
tially the same final temperature after cooling. We follow 
Bertschinger(1986) and we find that the shocked matter 
cools down to T ~ 10 5 K at z ~ 5. 



3. Number counts of primordial galaxies 

In order to quantify the global effect of the formation 
of primordial galaxies on the CMB, we apply our for- 
malism to a synthetic population of galaxies with masses 
1O 9 M < M < 1O 12 M . 

We first assume that the galaxy number density traces, 
within a linear bias, the abundance of collapsed DM ha- 
los, as predicted by the Press-Schechter (PS) mass func- 

We use an initial power- 



Press fc Schechter(1974) 



tion 

law spectrum with an effective spectral index n = — 2 
on galaxy scales. We express the amplitude of primor- 
dial matter fluctuations in terms of the rms variance in 



Viana & Liddle(1996), which corresponds to a bias factor 



b w 1.67. For fl = 1, the set of parameters corresponds to 
the "standard" biased cold DM model, which does fit nei- 



ther small- and large-scale velocities Vittorio et al.(1986) 



nor COBE normalisation. However, we take it as a study 
case for the computations, our second model is the low 
density cosmological model with Q = 0.3. In our picture 
(no stars are formed yet), the spheroid is gaseous. Its mass 
is related to the mass of the DM halo via M sp h = \Mh a i . 
The mass and luminosity of the central BH and thus the 
predicted SZ distortions are therefore inferred from M sph 
(Mbh = e BH M sph , and e BH = 2. lO" 3 ; cf §2). 

To compute the kinetic SZ term of a population 
of proto-galaxies, we need an estimate of their pecu- 
liar velocities with respect to the reference frame. As 



suggested by numerical simulations Bahcall et al.(1994) 
|Moscardini et al.(1996) , , we assume that velocities fol- 
low a Gaussian distribution. The peculiar velocity of each 
proto-galaxy is drawn from a Gaussian which is com- 
pletely defined by its rms value a v . In the range of red- 
shifts we have adopted, the structures arc in the linear 
regime, so that a v {z) — ctq/(z), where the redshift de- 



pendence of the velocities is given by /(z)Peebles(1980) 
|Peebles(1993l as a function of the cosmological pa- 
rameters. In this equation, a is the present-day rms 
peculiar velocity. It is related to the mass variance on 
m ass scale M, a(M) = (1 19» )("+ 3 )/ 6 a s M-^ n + 3 ^ 6 



Mathiesen &; Evrard(1998) , where n is the index of 
the power spectrum. The rms velocity can thus be com- 
puted for each mass scale. 



4. Results and discussion 

The shock-heated gas within proto-galaxies interacts with 
the CMB photons through the SZ effect (thermal and ki- 
netic) . These interactions generate secondary temperature 
anisotropics and spectral distortions. We simulate maps of 
the secondary anisotropies generated by a population of 
seeded proto-galaxies formed between redshift 5 and 10. 
The maps have a resolution of about 0.2 arcseconds to re- 
solve the galaxies and contain 600 x 600 pixels. The num- 
ber of sources of mass M at redshift z is derived from the 
PS mass function. Their positions are drawn at random 
in the map. The y and ST/T profiles for, respectively, the 
thermal and the kinetic effects are directly derived from 
the integration of the gas profile n e (R) along the line of 
ht (Eqs. | and H) assuming spherical symmetry. 
Similarly to the case of galaxy clusters, we assume that 
in the early stages of formation the gas settles into a hy- 
drostatic equilibrium within the DM potential. A univer- 



sal density profile is motivated by Navarro et al.(1996) 



However, the gas profile may be softer than that of the 
DM, and moreover the existence of a central cusp is "un- 
observed" Kravtsov et al.(1998) , We thus conservatively 
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adopt the following parametrised profile for the gas distri- 
bution: 



n e (R) = no 



, . R 



(J) 



where no is the central density, a is left as a free param- 
eter describing the steepness of the profile, whereas R c 
is identified with a core radius as in galaxy clusters. On 
cluster scales, R c is typically 10 to 30 times smaller than 
the cluster virial radius R V i r - In our model we introduce 
the parameter p = R V i r /R c which we vary, similarly to 
clusters, between 10 and 30. The central density no can 
be derived from the gas mass of the proto-galaxy using 
the following equation: 

M mh I ^ ] = m„M / n e (R) 4ttR 2 dR, (8) 



where the virial radius of the structure is given by: 
(GM) 1 /3 i 



(3^# ) 2/3 1 + W 



(9) 



for a critical universe. It is fixed solely by the mass and 
the collapse redshift z co a- 

We will give the results for a flat model with no cosmo- 
logical constant (Qo = 1) and an open model (Qq = 0.3). 
Varying the cosmological parameters will vary the number 
of proto-galaxies along the line of sight as well as their pe- 
culiar velocities. It will also modify their physical proper- 
ties, i.e. the size and velocity of the shock and thus the gas 
temperature. The two cosmological models represent the 
upper and lower bounds between which all other cosmo- 
logical models involving a non-zero cosmological constant 
fall. 

4-1. Compton distortion 

The CMB photons, scattering off the electrons of the 
ionised hot gas, induce a spectral distortion whose am- 
plitude is given by Eq. |l|. The FIRAS experiment has 
measured the mean Compton parameter resulting from 
all the interactions undergone by the photons. The result 
is Vfiras — l-^ 10~ 5 Fixscn ct al. (1996)] , . This stringent 
observational limit incorporates the (negligible) contribu- 
tion of the rather cold intergalactic medium and that of all 
other extragalactic signals. Among these signals, there is 
the contribution of the hot ionised gas in galaxy clusters. 
The global distortion induced by i ntra-clustcr gas has be en 



computed |Dc Luca ct al.(1995)| , |Barbosa et al.(1996) 



and found to be of the order of a few 10 -6 . In addition 
to galaxy clusters, one has to take into account the con- 
tribution of the proto-galaxy population in terms of the 
overall Compton distortion, y PG , induced by the scatter- 
ing of CMB photons on the shock- heated gas. 

Based on simulated maps, we predict y P Q and we com- 
pare it to the limit set by COBE-FIRAS. Among all the 



parameters of the model, there are four major quanti- 
ties that substantially affect the predictions of the mean 
Compton parameter. Two of them, a and p, are related 
to the gas distribution (Eq. Q). The two others are the 
fraction / of BH-seeded proto-galaxies and the BH-to- 
spheroid mass ratio cbh- We compare our predicted over- 
all distortion to the COBE-FIRAS limit and look for the 
combinations of parameters for which our predictions fit 
the observations. This allows us to constrain the assump- 
tions of our model. 

For a = 1/2, / = 100%, e BH = 2. 10~ 3 and both 
cosmological models, we find y P Q — 10~ 4 which ex- 
ceeds the observational value. In order for our prediction 
to be reconciled with the COBE-FIRAS limit, / must 
be only a few percent. This constraint on / strongly 



violates the actual observations Magorrian et al.(1998) 
|Richstone ct al.(1998)[ , . a = 1/2 is thus excluded by the 
limit on the global distortion whatever value we choose 
for p = R mr /R c . 

For a — 1 (i.e. an isothermal profile), an £7o = 1 
model, / = 100% and e B H = 2. 10 -3 , we find y PG > 
V firas whatever we adopt for p. The fraction / must 
be smaller than 75% for the prediction to be com- 
patible with the observational limit. Again, this frac- 
tion is significantly smaller than the 95% advocated by 



Magorrian et al.(1998). Such a constraint could rule out 



the isothermal profile. However, up to now, €bh was as- 
sumed to be constant and equal to 2. 10~ 3 . If we now 
use the lower limit of Magorrian et al.(1998), that is 
€bh = 10~ 3 , together with / = 95% or higher there is 
only a marginal agreement for p ~ 30 between the pre- 
dicted and measured distortions. In the open model case 
(f^o = 0.3), we find approximately the same results. For 
]Jpq to be compatible with y firas ^ P = 10, the fraction 
of BH seeded proto-galaxies should be smaller than 80% if 
e BH = 2. 10~ 3 or / < 95% if e BH = 1CT 3 . The predictions 
for p — 30 agree with observations in all the cases. 

For a — 3/2 (i.e. the gas profile approximates a King 
profile) and for both cosmological models, we find y P Q of 
about 10~ 6 to a few 10 -6 , a prediction compatible with 
V firas- This result remains valid for all values of p be- 
tween 10 and 30 including the highest boundaries of / and 
fflif, respectively 100% and 2. 10~ 3 . 

4-2. Predicting the angular power spectrum 

We choose the set of parameters associated with the 
isothermal profile which agrees with the COBE-FIRAS 
limit: a = 1, p = 30, / = 95% and e B H = 2. 10~ 3 . Within 
this context, we predict the upper limit on the contribu- 
tion to secondary temperature anisotropies induced by the 
SZ, thermal and kinetic effects, of the proto-galactic gas. 
We express this contribution in terms of an angular power 
spectrum plotted in figure |l] together with the main other 
well-known secondary anisotropies. 

At very small scales (I ~ a few 10 5 ) corresponding to 
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galactic scales, the kinetic SZ contribution of the shock- 
heated gas (Fig. [l] thick solid line for Q = 1 and thick 
dashed line for f2 = 0.3) is very large. It is interesting to 
note the good agreement between our results and those ob- 



tained by Peebles & Juszkicwicz(1998) for the scattering 
of the CMB photons by the cloudy proto-galactic plasma. 
The power spectrum of the kinetic SZ anisotropics for the 
Oq = 0.3 model is significantly larger than the flo = 1 
model. This is mainly due to the higher number of sources 
per unit comoving volume in open models. In all other 
flat cosmological models involving a non-zero cosmologi- 
cal constant, the power spectrum will lie between the two 
curves. The expected power spectrum due to the thermal 
SZ effect is not plotted in this figure. It is more than one 
order of magnitude smaller than the kinetic effect con- 
tribution. This is due to the efficiency of brcmsstrahlung 
which lowers the temperature down to a few 10 5 K. 

We compare the contribution of the proto-galaxies due 
to their SZ kinetic effect to the major sources of sec- 
ondary temperature anisotropies. In each case, we choose 
the most extreme cases for the comparison with our up- 
per limit prediction. The power spectra displayed in figure 
|l| are taken from the literature. The dotted line repre- 
sents the upper limit of the contribution of the inhomoge- 



neous reionisation as computed by Aghanim ct al.(1996) 



for a quasar lifetime of 10 7 yrs. The dot-dashed line repre- 
sents the Rees-Sciama effect Rees & Sciama(1968) , taken 

1, fl h = 0.25). The dashed 



from Seljak(1996) 



line represents the galaxy cluster contribution due to ki- 



netic SZ effect from Aghanim et al.(1998) with a cut-off 
mass of 10 14 Mq. The triple-dot-dashed line represents 



the Vishniac-Ostriker effect Ostriker & Vishniac(1986) 



Vishniac(1987)| , computed by |Hu fe Whitc(1996) 
total reionisation occurring at z 



ith a 



10. Finally, the solid 
thin line represents the power spectrum of the primary 
CMB anisotropies in a standard CDM model computed 
using the CMBFAST code [Seljak fe Zaldarriaga(1996)| , . 
The primary CMB anisotropies dominate at all scales 
larger than the damping around 5 arcminutes. At interme- 
diate scales, several effects take place among which the in- 
homogeneous reionisation, the Ostriker- Vishniac and the 
SZ effect. In figure [l], we do not plot the power spectra of 
the thermal SZ effect of galaxy clusters. It is about one 
order of magnitude larger than the kinetic SZ effect. At 
very small scales, the anisotropies are totally dominated 
by the proto-galactic contribution. 

5. Conclusions 

Previous studies on galaxy formation have evaluated the 
global distortion of a population of galaxies in the viri- 
alised regime. In these studies, the global Compton pa- 
rameter was found to be very small, and smaller than 
the COBE-FIRAS value. In contrast, our model focuses 
on a regime in which proto-galaxies undergo a BH forma- 
tion phase. During this phase, the proto-galactic matter is 




10 6 10 7 



Multipole 



Fig. 1. Power spectrum of the temperature anisotropies. The 
solid thin line represents the CMB primary anisotropies. The 
solid thick line represents the SZ kinetic effect of the proto- 
galaxies for Q = 1. The thick dashed line is obtained for 
Q = 0.3. The triple- dotted- dashed line represents the Vishniac- 
Ostriker effect. The dashed and dotted-dashed lines represent 
respectively the galaxy cluster contribution due to the kinetic 
SZ effect and the Rees-Sciama effect. 

shock-heated up to a few 10 8 K and cools down to 10 5 K. 
CMB photons undergo inverse Compton scattering on the 
heated gas. In addition, galaxy peculiar motions induce 
temperature anisotropies through the SZ kinetic effect. We 
have estimated the global Compton parameter due to a 
population of proto-galaxies and the expected power spec- 
trum of the induced secondary anisotropies. We find that 
there are four main parameters that control our model: the 
fraction / of BH-seeded proto-galaxies, the fraction cbh of 
the spheroid mass in the BH, the steepness of the density 
profile a and the gas core radius p = R V i r /R c . The com- 
parison between our predictions and the COBE-FIRAS 
observation constrains these parameters. Given the ob- 
served fraction of seeded galaxies, / = 95%, our results 
put rather strong constraints on the density profile and 
on cbh- Indeed, our predictions agree with the observa- 
tions whatever p if the density profile is an approximation 
to a King profile. On the contrary, if the density profile is 
isothermal, then the core radius must be at least 30 times 
smaller than the virial radius and the BH-to-spheroid mass 
ratio has to be small, of the order of 10 -3 . The computa- 
tions in the two extreme cosmological models show that 
the global Compton parameter due to proto-galaxies is 
not very sensitive to fio- 

We compare the power spectra of the different contri- 
butions to the temperature anisotropies. Our results show 
that the SZ effect of the very early shock-heated proto- 
galaxies could constitute the major source of CMB distor- 
tions on very small scales (arcsecond and sub-arcsecond 
scales). The anisotropies are likely to be detected and 
measured by future long baseline interferometers such as 
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ALMA. The shock heating is likely to contribute to the 
re-heating of the proto-galactic gas, which plays a role 



in galaxy formation theory. Blanchard et al.(1992) used 
preheating to modify the galaxy luminosity function, sup- 
pressing and finally delaying dwarf galaxy formation. We 
do not take into account this effect in our model, therefore, 
our results should be taken as an upper limit to the proto- 
galaxy contribution in terms of secondary anisotropies. 
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